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Abstract 

We present a determination of the pion-nucleon sigma-term based on a novel analysis of the 
ttN scattering amplitude in Lorentz covariant baryon chiral perturbation theory renormalized in 



<*S the extended-on-mass-shell scheme. This amplitude, valid up-to next-to-next-leading order in the 



chiral expansion, systematically includes the effects of the A(1232), giving a reliable description of 



the phase shifts of different partial wave analyses up to energies just below the resonance region. We 



obtain predictions on some observables that are within experimental bounds and phenomenological 
expectations. In particular, we use the center-of-mass energy dependence of the amplitude adjusted 
with the data above threshold to extract accurately the value of ct^n- Our study indicates that 
the inclusion of modern meson-factory and pionic-atom data favors relatively large values of the 
sigma term. We report the value a n N = 59(7) MeV. 

ON 

The sigma terms, a^N and a s , are observables of fundamental importance that embody the internal 
scalar structure of the nucleon, becoming an essential piece to understand the origin of the mass of 
the ordinary matter in the context of non-perturbative QCD. The pion-nucleon sigma term, a n N, also 
plays a role in the study of the equation of state of the nuclear and neutronic systems pQ and is 
a key ingredient in investigations of the QCD phase diagram that explore the restoration of chiral 
^ ■ symmetry in cold and dense matter [2]. On the other hand, a n N and a s , as properties quantifying the 
^ ■ response of the nucleons as probed by scalar interactions, appear in the hadronic matrix elements of 
the spin-independent neutralino- nucleon elastic scattering cross section [31 HI [5]. Unfortunately, our 
current knowledge of these quantities is far from satisfactory and they have become the main source of 
uncertainty in the interpretation of experimental searches of supersymmetric dark matter [31 H]. With 
the advent of experimental results on direct searches of dark matter, Ellis et al recently pled for a more 
accurate experimental determination of the sigma terms [4]. This has triggered an intense campaign 
for the obtention of these matrix elements from first principles using LQCD simulations [5J, although a 
revision of the experimental discrepancies is still missing. In this letter, we focus on the extraction of 
the a n N from nN scattering data. 

A connection between the ttN scattering amplitudes and a n N can be established, using the chiral 
symmetry of the strong interactions, at the Cheng-Dashen point [7J [8], which lies in the unphysical 
region of the amplitude (W c d = y/s^d = Mn, t C( i = 2m\ with the physical region above W t h = m n + 
and t < 0). At this point, one has E = a n N + + Ar, where £ is proportional to the Born- 
subtracted isoscalar nN scattering amplitude and the pion semileptonic decay constant squared f%, 
A CT ~ 15 MeV [9] and Ar is the remainder of the relation which chiral symmetry constrains to be 
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Table 1: Main physical observables obtained from the N 2 LO nN scattering amplitude in the EOMS 
renormalization scheme fitted to different PW analyses up to W max = 1.2 GeV (W max = 1.16 GeV for 
EM). Results of the scattering lengths are in units of 10~ 2 m" 1 . 

Xd.o.f, h A g nN A GT [%} qp + 

KH [EE] 0.75 3.02(4) 13.51(10) 4.9(8) -1.2(8) 8.7(2) 
GW [12J 0.23 2.87(4) 13.15(10) 2.1(8) -0.4(7) 8.2(2) 
EM [22] 0.11 2.99(2) 13.12(5) 1.9(4) 0.2(3) 7.7(1) 



~ 0(m\) and negligible (A# < 1 MeV) JTU]. The conventional way to perform the extrapolation 
from the physical to the subthreshold region is by means of an energy-dependent parameterization of 
the experimental data in partial waves (PW) supplemented by dispersion relations that impose strong 
analyticity and unitarity constraints from the high-energy data onto the low-energy scattering amplitude 
[7J. There are two of such analyses to be considered in this letter: those of the Karlsruhe- Helsinki 
(KH) [TTJ and the George- Washington [T2] (GW) groups. In particular, KH reports a £ = 64(8) MeV 
that leads to a n N = 49(8) MeV [IT], while the analysis of Gasser et al, based on the KH amplitudes, 
gives the classical result a^N — 45 MeV [9j. On the other hand, the analyses of the GW group, including 
modern pion factory data, obtain larger values, S = 79(7) MeV that leads to ovjv = 64(7) MeV [T3] . 

The 7r N scattering observables can be studied using chiral perturbation theory (%PT), which is 
the effective field theory of QCD at low energies [HJ US]. The advantage of %PT over the dispersive 
methods discussed above stems from the fact that this approach imposes the chiral Ward identities of 
QCD valid up to a certain order in the chiral power counting. Besides providing the corrections to the 
relation at the Cheng-Dashen point, it establishes a more direct connection between the isoscalar nN 
scattering amplitude and cr^-N- I n feet, up to 0(p 3 ) both observables are connected by the same low- 
energy constant (LEC) c\ [TBI HZ]- The LECs are parameters that accompany the effective operators 
in the chiral Lagrangian and are fixed by matching to QCD or by comparison with experimental data. 
Therefore, in one can predict a n N fixing the relevant LEC using nN scattering data and avoiding 
the extrapolation to the Cheng-Dashen point. 

In this note, we report on the extraction of a n N from irN scattering data and using a chiral represen- 
tation scattering amplitude that is obtained in Lorentz covariant B^PT renormalized in the extended- 
on-mass-shell (EOMS) scheme. A detailed explanation of the advantages of this approach as compared 
with previous works has been presented by J.M. Alarcon [T5] in his talk at this School. We have per- 
formed a calculation up to N 2 LO accuracy explicitly including the effects of the A (1232) resonance in 
the 5-counting, that at low-energies and below the resonance region, takes into account the fact that the 
diagrams with the A are suppressed in comparison with those with the nucleon |19] . A more thorough 
presentation of our results can be found in [20J . In the following, we outline some results of the resulting 
chiral representation on main observables and more specificaly on the pion-nucleon sigma term. 

The calculation in the EOMS scheme [21] proceeds as explained in [18]. On top of that, we include 
Born-terms with an intermediate A(1232), whereas the corresponding loops with A propagators are of 
higher-order. We then fix the values of the LECs fitting the CM energy dependence of the 2 S- and 
4 P-wave phase shifts obtained from the chiral amplitude to the latest solutions of the KH [TTJ and 
GW [12] groups. In addition, we include the analysis of the Matsinos' group (EM) [22] which focuses 
in the PW parameterization of the data at very low energies without imposing dispersive constraints 
from the high-energy region. The fits are done from the lowest CM energies above threshold, Wth — 
1.08 GeV, up to W max = 1.2 GeV which is below the A(1232) region (the EM analysis only reaches 
W max ~ 1.16 GeV). The parameters g A , M N , Ma, m n and f n are fixed to their experimental values 
g A = 1.267, M N = 939 MeV, M A =1232 MeV, m w = 139 MeV and = 92.4 MeV. We work with 
physical values which is equivalent to a reordering of the chiral series [23]. The iVA axial coupling can 
be determined using the A(1232) width, Ta = 118(2) MeV, giving h A = 2.90(2) [19]. However, we also 
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Figure 1: (Color on-line) Phase shifts given by the Lorentz covariant N 2 LO nN scattering amplitude 
in the EOMS scheme fitted to the GW solution (circles) [12] up to W max = 1.2 GeV. 

fit Ha to compare the result with this value. 

The main results on physical observables are shown in Table [TJ The errors quoted there are only of 
statistical origin and additional theoretical uncertainties are to be added. We plot in Fig. [1] the phase 
shifts of the S- and P-waves given by the N 2 LO ttN scattering amplitude in the EOMS scheme fitted to 
the GW solution (circles) p2] up to W max = 1.2 GeV. Similar plots can be obtained for the KH and EM 
solutions. The figure shows that the description for the lowest PWs is very good up to energies below 
the A(1232) resonance. This is reflected by the small x^of listed in the second column of Table [TJ 
showing that the description of GW and EM PW analyses is better. It is also worth remarking that the 
problems encountered by the IR prescription to describe the P\\ wave of the GW group [23] disappear 
in the present analysis. 

As one can infer from the third column of Table [TJ only the GW solution gives a result on 
that is perfectly compatible with the determination from the A(1232) width. In the fourth column, 
we show the values obtained for the nN coupling that, compared with the axial coupling g^, gives 
the GT discrepancy Act in the fifth column (see Ref. [23] for details). One can see that the results 
extracted from the KH, GW and EM PWs are completely consistent with the values reported by 
these collaborations, g vN = 13.41(14) [llj, g nN = 13.15(1) [12J and g vN = 12.98(12) [22] respectively. 
However, the KH result is in clear disagreement with the numbers independently extracted from NN- 
scattering (g n N — 13.0) [21] and pionic atom (g^N = 13.12(9)) [25] data, leading to a Act that may 
be considered too large p2]. Finally, our full results, with the explicit inclusion of the A, confirm the 
findings shown in [18]. Namely, that the large violation of the GT discrepancy obtained in IR (~20%) 
is a spurious effect introduced by this scheme rather than a problem of chiral convergence in the irN 
system. 

Results for the isoscalar (ag + ) and isovector (% + ) scattering lengths are shown in the last two 
columns of Table [TJ These scattering parameters are interesting because they can be compared with 
the values independently extracted from pionic-atom data. Moreover, parameter enclosing 
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Table 2: Values of the 0(p 2 ) LECs in units of GeV -1 and of a^N i n MeV obtained from the different 
ttN PW analyses. 

£l £2 £3 £4 OV/y 

KH -0.80(6) 1.12(13) -2.96(15) 2.00(7) 43(5) 
GW -1.00(4) 1.01(4) -3.04(2) 2.02(1) 59(4) 
EM -1.00(1) 0.58(3) -2.51(4) 1.77(2) 59(2) 



information on the isoscalar scattering amplitude, is tightly connected with a^N 0- I n this sense, the 
latest analyses of pionic-atom data clearly favor positive values for this observable, = 0.0076(31) 
m^ 1 |25j. As we can see by direct comparison with the values shown in Table [IJ only the GW and 
EM analyses are compatible with ajt. > 0. The impact that the pionic-atom result for has on the 
value of a nN is addressed below. In case of the isovector length, KH and GW are consistent with the 
pionic-atom value Oq + = 0.0861(9) [23]. For a thoughtful discussion on the discrepancy of the isovector 
scattering length extracted from the low-energy scattering data by the EM analyses, = 0.0797(11) 
m" 1 , see Refs. [22]. 

We have explored the subthreshold region and obtained the subthreshold coefficients which result 
from a Taylor expansion of the amplitudes about the point W c d and t = [7]. In general, we observe 
a good agreement on the leading coefficients as compared with those reported by the KH and GW 
solutions, what positively speaks of the potential of our amplitude to connect, in the CM energy 
W, the physical region to the Cheng-Dashen point. However, it is known that, at N 2 LO, the chiral 
representation underestimates the extrapolation to t — 2m 2 which is dominated by the threshold to two 
pions [HI ELY] • Nonetheless, and as was mentioned above, B%PT can predict a n N exclusively from the 
dependence in the CM energy W via the LEC C\ fitted to scattering data and avoiding the extrapolation 
into the unphysical region. 

The expression for a n N can be obtained either from the scalar form factor or the quark-mass depen- 
dence of the nucleon mass. The result in EOMS-BxPT is 



o 3g\ m l ( 3 M n - m l . m, 



(TnN = -4cim_ - — oro * , = = arccos — — + m n log 



16n 2 f 2 M N \^4M 2 N -ml 2M N & 

Eq. ([1]) is valid up to NLO accuracy in the chiral expansion and has an uncertainty from higher-order 
contributions, which we estimate computing the next subleading correction coming, in the ^-counting, 
from the loop with an insertion of the A [19]. This amounts to a contribution of ~ —6 MeV to be 
compared with the NLO term of —19 MeV. It is important to stress that we take this correction as 
an irreducible uncertainty in our determination. In order to explicitly add this contribution one has to 
include the same type of terms, arising from loops with one insertion of a A propagator, into the nN 
scattering amplitude studying the changes produced in the LECs. 

To obtain the values of ovat, we consider fits to the KH and GW PW phase shifts and for various 
W max , from 1.14 GeV to 1.2 GeV in intervals of 0.01 GeV. For the analysis of EM we perform the 
same study up to its maximum CM energy ~ 1.16 GeV. The purpose of this strategy is to take into 
account the dispersion of ct^n against the data set included in the fits. In Table |2l we show the mean 
and mean deviation for c\ and a n N- We also include results for the other 0(p 2 ) LECs for comparison 
purposes. In general, we find that for the KH and GW analyses ovjv decreases from 48 MeV to 39 MeV 
and from 65 MeV to 54 MeV respectively, when we increase W max within the given bounds. The 
results for the EM solution remain quite stable. As we can see, the values of a W N are not completely 
consistent among each other. The number obtained from the KH PWs agrees with the canonical result 
a n N — 45 MeV, whereas those of the GW and EM solutions completely agree with each other, and 
with the one determined by the GW group, in a value ~ 16 MeV larger. 
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The consistency between the results derived from the GW and EM solutions is very remarkable 
since these are quite different analyses having both in common the implementation of the wealth of 
low-energy data collected along the last 20 years in meson factories [121 122] with many points not 
included in KH [llj. Besides that, the positive values determined from pionic-atom data for aj, , which 
are consistent with those extracted from the GW and EM PWs, strongly constrain a n ^, raising it for 
>7 MeV as compared with the KH values [131 M- We also stress that only our results based in the GW 
PW analysis are perfectly compatible with all the phenomenology discussed. With these considerations, 
one may obtain the following value for a^N, as it is extracted from the analysis of nN modern scattering 
data and using B%PT, 

a vN = 59(7)MeV. (2) 

The error includes the higher-order uncertainties estimated above added in quadrature with those given 
by the dispersion of the values in the average of the GW and EM results. If one were to include the 
KH result in this estimation, the result would by slightly reduced by 2-3 MeV. 

As a concluding observation we want to address the fact that this large result on a^N might be 
in conflict with some established phenomenology. In particular, it resurrects an old puzzle concerning 
the strangeness content of the nucleon [13]. This relies on the relation that is obtained in HB%PT up 
to NLO accuracy among the S , [/(3)i?-breaking of the baryon-octet masses, ct^n and the observable y 
quantifying the strangeness content of the nucleon [2S] • For the value of the sigma term obtained in the 
present work, this relation leads to a contribution of the strange quark to the nucleon mass of several 
hundreds of MeV. This is a physical scenario that is very hard to understand. In this regard, it has 
been recently shown that the HB framework does not provide a reliable description of iS[/(3)i?-breaking 
phenomenology and, in particular, of the quark- mass dependence of the baryon octet masses [27]. In 
fact, in the analysis of LQCD data done using Lorentz covariant B^PT (EOMS scheme) and including 
the decuplet-baryon contributions, it has been obtained that a a n N of ~60 MeV is perfectly compatible 
with a negligible strangeness in the nucleon [27j . Therefore, and in our opinion, the relation between 
the 5 , [/(3)i?-splitting of the baryon octet masses and the sigma terms ought to be revised. 

Another caveat for a relatively large value of the pion-nucleon sigma term arises in chiral approaches 
of nuclear matter [Tj |2] , in which pa n N controls the leading contribution in the density (p) dependence 
of the quark condensate. At this order, a large value for a^N — 60 MeV would imply a vanishing in- 
medium quark condensate at ~ 2p Q (with p the nuclear saturation density), while this occurs at ~ 3po 
using the canonical value of 45 MeV [TJ. However, a non- vanishing value of the in- medium temporal 
component of the pion axial coupling, f t , is also a necessary condition for the spontaneous chiral 
symmetry breaking [2H]- The leading in- medium contribution to f t is controlled by the combination 
c 2 + c 3 , which is not directly related to ct^tv at this order. The difference in c 2 + c 3 between KH and 
GW sets in Table [2] is only of 10%, which means that the running with density for f t , as obtained 
in [2H] (which also includes iViV interactions), differs also on a 10%. This figure is much smaller than 
the difference of around a 30% in the density dependence of the quark condensate due to the different 
values of a n N [21], so that the vanishing of f t occurs at ~ S.5p for both PW analyses. A more thorough 
analysis is then necessary to put in agreement the simultaneous vanishing of f t and the quark condensate 
with density in order to properly discuss about chiral symmetry restoration in nuclear matter. 

In summary, we have presented the values on a n N resulting from a novel analysis of the nN scattering 
amplitude in Lorentz covariant B%PT within the EOMS scheme up to 0(p 3 ) [18] and including the effects 
of the A(1232) explicitly in the (^-counting. We have found that we perfectly describe the PW phase 
shifts of the KH, GW and EM groups up to energies below the A-resonance region. Our amplitudes 
are suitable to extract an accurate value of a n ^ from scattering data and avoiding the extrapolation to 
the unphysical region. Our work confirms the results obtained from dispersion relations by the KH and 
GW groups. This agreement, in turn, gives us further confidence in our analysis based on B%PT. We 
then ratify the discrepancy between them and give further support to the GW result after studying the 
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latest EM solution. We conclude that modern data lead to a relatively high value of cr^jv- 
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